cd "/Users/mmoral/Dropbox/Sabancı University/-ONGOING/On the Relationship between Party Polarization and Citizen Polarization/"
clear all
use "Replication Data.dta"
xtset iso elec

** Figures 1 and 2
*set scheme mmoral3															/* Scheme file is available from the authors upon request */
twoway 	(lfit resleftvot date, lcolor(black) lpattern(dot) yaxis(1)) (lfit resleftabs date, lcolor(gs6) lpattern(shortdash) yaxis(1)) ///
		(scatter resleftvot date, msymbol(smcircle) mcolor(black) yaxis(1)) (scatter resleftabs date, msymbol(smcircle_hollow) mcolor(gs6) yaxis(1)) ///
		(lowess uwedis date, lcolor(black) lpattern(solid) lwidth(vthin) msymbol(i) yaxis(2)), ///
		by(iso, note("")) legend(order(3 "Voters (Observed)" 1 "Voters (Linear Fit)" 5 "Party Polarization (Lowess)" 4 "Non-voters (Observed)" 2 "Non-voters (Linear Fit)") rows(2) symxsize(3) size(vsmall)) ytitle("Citizen Polarization", axis(1)) /// 
		xtitle("Election Year") xlab(1970(10)2020) xsca(r(1968 2020)) ylab(1(.5)3.5, axis(1)) ylab(1(.5)3.5, axis(2)) name(vote_all, replace) 
		
twoway 	(lfit resleftsop1 date, lcolor(black) lpattern(dot) yaxis(1)) (lfit resleftsop2 date, lcolor(gs6) lpattern(shortdash) yaxis(1)) (lfit resleftsop3 date, lcolor(gs12) lpattern(longdash) yaxis(1)) ///
		(scatter resleftsop1 date, msymbol(smcircle) mcolor(black) yaxis(1)) (scatter resleftsop2 date, msymbol(smcircle_hollow) mcolor(gs6) yaxis(1)) (scatter resleftsop3 date, msymbol(smsquare) msize(*.75) mcolor(gs12) yaxis(1)) ///
		(lowess uwedis date, lcolor(black) lwidth(vthin) msymbol(i) yaxis(2)), ///
		by(iso, note("")) legend(order(4 "High Sophistication (Observed)" 1 "High Sophistication (Linear Fit)" 5 "Med. Sophistication (Observed)" 2 "Med. Sophistication (Linear Fit)" 6 "Low Sophistication (Observed)" 3 "Low Sophistication (Linear Fit)" 7 "Party Polarization (Lowess)") rows(2) symxsize(3) size(vsmall)) ytitle("Citizen Polarization", axis(1)) /// 
		xtitle("Election Year") xlab(1970(10)2020) xsca(r(1968 2020)) ylab(1(.5)3.5, axis(1)) ylab(1(.5)3.5, axis(2)) name(soph_all, replace)
		
preserve 
keep if inlist(iso, 1036, 1208, 1276, 1528, 1752, 1840)
twoway 	(lfit resleftvot date, lcolor(black) lpattern(dot) yaxis(1)) (lfit resleftabs date, lcolor(gs6) lpattern(shortdash) yaxis(1)) ///
		(scatter resleftvot date, msymbol(smcircle) mcolor(black) yaxis(1)) (scatter resleftabs date, msymbol(smcircle_hollow) mcolor(gs6) yaxis(1)) ///
		(lowess uwedis date, lcolor(black) lpattern(solid) lwidth(vthin) msymbol(i) yaxis(2)), ///
		by(iso, rows(2) note("")) legend(order(3 "Voters (Observed)" 1 "Voters (Linear Fit)" 5 "Party Polarization (Lowess)" 4 "Non-voters (Observed)" 2 "Non-voters (Linear Fit)") rows(2) symxsize(3) size(small)) ytitle("Citizen Polarization", axis(1)) /// 
		xtitle("Election Year") xlab(1970(10)2020) xsca(r(1968 2020)) ylab(1(.5)3.5, axis(1)) ylab(1(.5)3.5, axis(2)) name(vote_ind, replace)
twoway 	(lfit resleftsop1 date, lcolor(black) lpattern(dot) yaxis(1)) (lfit resleftsop2 date, lcolor(gs6) lpattern(shortdash) yaxis(1)) (lfit resleftsop3 date, lcolor(gs12) lpattern(longdash) yaxis(1)) ///
		(scatter resleftsop1 date, msymbol(smcircle) mcolor(black) yaxis(1)) (scatter resleftsop2 date, msymbol(smcircle_hollow) mcolor(gs6) yaxis(1)) (scatter resleftsop3 date, msymbol(smsquare) msize(*.75) mcolor(gs12) yaxis(1)) ///
		(lowess uwedis date, lcolor(black) lwidth(vthin) msymbol(i) yaxis(2)), ///
		by(iso, rows(2) note("")) legend(order(4 "High Sophistication (Observed)" 1 "High Sophistication (Linear Fit)" 5 "Med. Sophistication (Observed)" 2 "Med. Sophistication (Linear Fit)" 6 "Low Sophistication (Observed)" 3 "Low Sophistication (Linear Fit)" 7 "Party Polarization (Lowess)") rows(2) symxsize(3) size(vsmall)) ytitle("Citizen Polarization", axis(1)) /// 
		xtitle("Election Year") xlab(1970(10)2020) xsca(r(1968 2020)) ylab(1(.5)3.5, axis(1)) ylab(1(.5)3.5, axis(2)) name(soph_ind, replace) nodraw
restore

loc i=1
foreach graph in vote_ind soph_ind vote_all soph_all{
graph disp `graph', scale(.9)
gr_edit .plotregion1.subtitle[1].style.editstyle size(medium) editcopy 
if `i'<=2 {
graph export "Graphs/Figure`i'.pdf", replace
}
if `i'>2 {
loc j=`i'-2
graph export "Graphs/App_Figure`j'.pdf", replace
}
loc i=`i'+1
}

** Analyses
clear all
use "Replication Data.dta"
xtset iso elec

** Table 1
eststo m1: sureg (d.uwedis l.uwedis d.resleft l.resleft) (d.resleft l.resleft d.uwedis l.uwedis), dfk corr
eststo m2: sureg (d.uwedis l.uwedis d.resleftvot l.resleftvot) (d.uwedis l.uwedis d.resleftabs l.resleftabs), dfk corr
eststo m3: sureg (d.uwedis l.uwedis d.resleftsop1 l.resleftsop1) (d.uwedis l.uwedis d.resleftsop2 l.resleftsop2) (d.uwedis l.uwedis d.resleftsop3 l.resleftsop3), dfk corr
eststo m4: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis) (d.resleftabs l.resleftabs d.uwedis l.uwedis), dfk corr
eststo m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis), dfk corr

* Long-run multipliers
forval i=1/5{
est restore m`i'
mat b=e(b)
mat V=e(V)
foreach j of numlist 1 5 9{
sca x`j'=b[1,`=`j'+2']
sca y`j'=b[1,`j']
sca varx`j'=V[`=`j'+2',`=`j'+2']
sca vary`j'=V[`j',`j']
sca covx`j'y`j'=V[`=`j'+2',`j']
sca b`j'_LRM=x`j'/-y`j'
sca se`j'_LRM=sqrt((1/y`j'^2)*varx`j'+(x`j'^2/y`j'^4)*vary`j'-2*(x`j'/y`j'^3)*covx`j'y`j')
sca t`j'_LRM=b`j'_LRM/se`j'_LRM
estadd sca b`j'_LRM
estadd sca se`j'_LRM
disp "p`j'_LRM=" %9.3f ttail(e(dfk2_adj), t`j'_LRM)
}
}
/*
esttab m1 m2 m3 m4 m5 using "Graphs/Table1.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b5_LRM b9_LRM se1_LRM se5_LRM se9_LRM, layout(@ @ @ "@ @ @" "@ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) ///
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.1" "Model.2" "Model.3" "Model.4" "Model.5") order(L.uwedis D.resleft L.resleft D.uwedis) ///
rename(D.resleftvot D.resleft D.resleftabs D.resleft D.resleftsop1 D.resleft D.resleftsop2 D.resleft D.resleftsop3 D.resleft L.resleftvot L.resleft L.resleftabs L.resleft L.resleftsop1 L.resleft L.resleftsop2 L.resleft L.resleftsop3 L.resleft) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" _cons Constant) substitute(\_ _) fragment
*/

** In-text descriptive stats --i.e., total number of parties in examined elections
egen totpar=total(pidN)
sum totpar
disp `r(mean)'/`r(N)'
disp `r(N)'/19

** In-text explanation of the substantive significance of the findings --i.e., effect size
est restore m4
test [D_resleftvot]L.uwedis = [D_resleftabs]L.uwedis
mat lis e(V)

est restore m5
test [D_resleftsop1]L.uwedis
test [D_resleftsop1]L.uwedis = [D_resleftsop3]L.uwedis
test [D_resleftsop1]D.uwedis = [D_resleftsop3]D.uwedis

disp sqrt(.00216644+.00537254-2*.00157654)
disp (([D_resleftsop1]D.uwedis-[D_resleftsop3]D.uwedis)/.06622613)^2

** In-text explanation of the difference of estimated coefficients
est restore m1
mat b=e(b)
sum D.resleft if e(sample)
disp `r(sd)'*_b[D_uwedis:D.resleft]
sum L.resleft if e(sample)
disp `r(sd)'*_b[D_uwedis:L.resleft]
sum D.uwedis if e(sample)

sum L.uwedis if e(sample)
disp `r(sd)'*_b[D_resleft:L.uwedis]
sum D.uwedis if e(sample)
disp `r(sd)'*_b[D_resleft:D.uwedis]
sum D.resleft if e(sample)

est restore m5
estat sum
sum D.uwedis if e(sample)
disp `r(sd)'*_b[D_resleftsop1:D.uwedis]
sum D.resleftsop1 if e(sample)

** In-text descriptive stats --i.e., average time between consecutive election cycles
gen year_dif=year-L.year if iso==L.iso
tabstat year_dif if e(sample), by(countryname) stats(mean N)
disp 3.832258 *5
tabstat date, by(countryname) stats(min max N)

** Figure 3
preserve
clear all
use "Replication Data.dta"
xtset iso elec

eststo m1: sureg (d.uwedis l.uwedis d.resleft l.resleft) (d.resleft l.resleft d.uwedis l.uwedis), dfk corr 																	
eststo m2: sureg (d.resleft l.resleft d.uwedis l.uwedis) (d.uwedis l.uwedis d.resleft l.resleft), dfk corr 																	/* Equations are re-ordered for the loop below, Model.1 */
eststo m3: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis) (d.resleftabs l.resleftabs d.uwedis l.uwedis), dfk corr
eststo m4: sureg (d.resleftabs l.resleftabs d.uwedis l.uwedis) (d.resleftvot l.resleftvot d.uwedis l.uwedis), dfk corr														/* Equations are re-ordered for the loop below, Model.4 */
eststo m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis), dfk corr
eststo m6: sureg (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis) (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis), dfk corr	/* Equations are re-ordered for the loop below, Model.5 */

forval j=1/6{
est restore m`j'
mat b`j'=e(b)
forval i=1/3{
sca b`j'`i'=b`j'[1,`i']
}
sca lrm`j'=b`j'3/-b`j'1
sca t`j'1=b`j'2
sca t`j'2=(1+b`j'1)*b`j'2+(b`j'3-b`j'2)
forval i=3/6{
sca t`j'`i'=(1+b`j'1)*t`j'`=`i'-1'
}
forval i=1/6{
sca eff`j'`i'=t`j'`i'/lrm`j'
}
sca toteff`j'1=eff`j'1
forval i=2/6{
sca toteff`j'`i'=toteff`j'`=`i'-1'+eff`j'`i'
}
}

forval j=1/6{
gen lrm`j'=lrm`j'
forval i=1/6{
gen t`j'`i'=t`j'`i'
gen eff`j'`i'=eff`j'`i'
gen toteff`j'`i'=toteff`j'`i'
}
}

keep lrm1-toteff66
keep in 1
gen id=1

rename lrm* lrm_m*
foreach var in eff toteff t{
forval i=1/6{
rename (`var'`i'1 `var'`i'2 `var'`i'3 `var'`i'4 `var'`i'5 `var'`i'6) (`var'1_m`i' `var'2_m`i' `var'3_m`i' `var'4_m`i' `var'5_m`i' `var'6_m`i')
}
}

reshape long lrm_m t1_m t2_m t3_m t4_m t5_m t6_m eff1_m eff2_m eff3_m eff4_m eff5_m eff6_m toteff1_m toteff2_m toteff3_m toteff4_m toteff5_m toteff6_m, i(id) j(model)
rename *_m *
reshape long t eff toteff, i(model) j(period)

lab def models 1 `""Model.1 - {&Delta}Party Polarization" " ""' 2 `""Model.1 - {&Delta}Citizen Polarization" "(All Respondents)""' 3 `""Model.4 - {&Delta}Citizen Polarization" "(Voters)""' 4 `""Model.4 - {&Delta}Citizen Polarization" "(Non-voters)""' 5 `""Model.5 - {&Delta}Citizen Polarization" "(High Sophistication)""' 6 `""Model.5 - {&Delta}Citizen Polarization" "(Med. Sophistication)""', replace
lab val model models
replace eff=.a if period==1 & !inlist(model, 1, 2, 5)	/* Short-term effects of citizen polarization and party polarization are insignificant in other models and equations */
twoway (bar eff period, barwidth(.65) fcolor(black%35) lcolor(black%35)) (bar eff period if period==1 & inlist(model, 1, 2, 5), barwidth(.65) fcolor(black%65) lcolor(black%65)), xtitle(" ") xlab(1 "t{subscript:0}" 2 "t{subscript:1}"  3 "t{subscript:2}" 4 "t{subscript:3}" 5 "t{subscript:4}" 6 "t{subscript:5}")  by(model, note("{it:Note:} Short-run, instantaneous effects are indicated with darker gray.") imargin(medsmall) /*hole(4)*/ legend(off) rows(2) /*hole(3)*/) subtitle(, size(small)) ///
ytitle("Short- and Long-term Effects", size(small)) ylab(0(.25).75, gmax)
graph export "Graphs/Figure3.pdf", replace
restore

** Appendices
** Appendix B - Descriptive statistics
tab decade, gen(d)
tab pr_rule, gen(pr)
qui sureg (d.uwedis l.uwedis d.resleft l.resleft) (d.resleft l.resleft d.uwedis l.uwedis), dfk corr
eststo desc1: estpost tabstat deluwedis delresleft delresleftvot delresleftabs delresleftsop1 delresleftsop2 delresleftsop3 ///
luwedis lresleft lresleftvot lresleftabs lresleftsop1 lresleftsop2 lresleftsop3 ///
delwedis delpolv delpi delestray delsartori lwedis lpolv lpi lestray lsartori delresleftpkn delresleftpuk delreslefteduch delreslefteducl lresleftpkn lresleftpuk lreslefteduch lreslefteducl d2 d3 d4 d5 delenep lenep delav_dismag lav_dismag pr2 pr3 mdif other_scale_d compulsory resleftsop_alt1 resleftsop_alt2 resleftsop_alt3 resleftsop_alt4 if e(sample), stats(min mean max sd N) columns(statistics)
preserve
drop if gap==1
qui sureg (d.uwedis l.uwedis d.resleft l.resleft) (d.resleft l.resleft d.uwedis l.uwedis), dfk corr
eststo desc2: estpost tabstat deluwedis delresleft delresleftvot delresleftabs delresleftsop1 delresleftsop2 delresleftsop3 ///
luwedis lresleft lresleftvot lresleftabs lresleftsop1 lresleftsop2 lresleftsop3 if e(sample), stats(min mean max sd N) columns(statistics)
restore
preserve
keep if pr_rule<2
qui sureg (d.uwedis l.uwedis d.resleft l.resleft) (d.resleft l.resleft d.uwedis l.uwedis), dfk corr
eststo desc3: estpost tabstat deluwedis delresleft delresleftvot delresleftabs delresleftsop1 delresleftsop2 delresleftsop3 ///
luwedis lresleft lresleftvot lresleftabs lresleftsop1 lresleftsop2 lresleftsop3 if e(sample), stats(min mean max sd N) columns(statistics)
restore
qui sureg (d.uwedis l.uwedis d.resleft l.resleft l2.uwedis) (d.resleft l.resleft d.uwedis l.uwedis l2.resleft), dfk corr
eststo desc4: estpost tabstat deluwedis delresleft delresleftvot delresleftabs delresleftsop1 delresleftsop2 delresleftsop3 ///
luwedis lresleft lresleftvot lresleftabs lresleftsop1 lresleftsop2 lresleftsop3 if e(sample), stats(min mean max sd N) columns(statistics)
/*
esttab desc1 using "Graphs/Appendix B.tex",tex replace cells("mean(fmt(3)) sd(fmt(3)) min(fmt(3)) max(fmt(3)) count(fmt(0))") collabels("Mean" "Std.Dev." "Min." "Max." "N" ) noobs nonumber label title("Descriptive Statistics") longtable fragment	
esttab desc2 using "Graphs/Appendix B.tex",tex append cells("mean(fmt(3)) sd(fmt(3)) min(fmt(3)) max(fmt(3)) count(fmt(0))") collabels("Mean" "Std.Dev." "Min." "Max." "N" ) noobs nonumber label title("Descriptive Statistics") longtable fragment
esttab desc3 using "Graphs/Appendix B.tex",tex append cells("mean(fmt(3)) sd(fmt(3)) min(fmt(3)) max(fmt(3)) count(fmt(0))") collabels("Mean" "Std.Dev." "Min." "Max." "N" ) noobs nonumber label title("Descriptive Statistics") longtable fragment
esttab desc4 using "Graphs/Appendix B.tex",tex append cells("mean(fmt(3)) sd(fmt(3)) min(fmt(3)) max(fmt(3)) count(fmt(0))") collabels("Mean" "Std.Dev." "Min." "Max." "N" ) noobs nonumber label title("Descriptive Statistics") longtable fragment
*/
  
** Robustness and Sensitivity Checks
** Appendix C - Other Polarization Measures
clear all
use "Replication Data.dta"
xtset iso elec

foreach parpol of varlist wedis polv pi estray sartori{
eststo m1_`parpol': qui sureg (d.`parpol' l.`parpol' d.resleft l.resleft) (d.resleft l.resleft d.`parpol' l.`parpol'), dfk corr
eststo m2_`parpol': qui sureg (d.`parpol' l.`parpol' d.resleftvot l.resleftvot) (d.`parpol' l.`parpol' d.resleftabs l.resleftabs), dfk corr
eststo m3_`parpol': qui sureg (d.`parpol' l.`parpol' d.resleftsop1 l.resleftsop1) (d.`parpol' l.`parpol' d.resleftsop2 l.resleftsop2) (d.`parpol' l.`parpol' d.resleftsop3 l.resleftsop3), dfk corr
eststo m4_`parpol': qui sureg (d.resleftvot l.resleftvot d.`parpol' l.`parpol') (d.resleftabs l.resleftabs d.`parpol' l.`parpol'), dfk corr
eststo m5_`parpol': sureg (d.resleftsop1 l.resleftsop1 d.`parpol' l.`parpol') (d.resleftsop2 l.resleftsop2 d.`parpol' l.`parpol') (d.resleftsop3 l.resleftsop3 d.`parpol' l.`parpol'), dfk corr
}

foreach parpol of varlist wedis polv pi estray sartori{
forval i=1/5{
est restore m`i'_`parpol'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 5 9{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
}

/*
local i=1
esttab m`i'_wedis m`i'_polv m`i'_pi m`i'_estray m`i'_sartori using "Graphs/Appendix C`i'.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM se1_LRM b5_LRM se5_LRM b9_LRM se9_LRM, fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress ///
rename(L.polv L.wedis L.pi L.wedis L.estray L.wedis L.sartori L.wedis D.polv D.wedis D.pi D.wedis D.estray D.wedis D.sartori D.wedis) ///
coeflabels(D.resleft "$\Delta\$Citizen Polarization" L.resleft "Citizen Polarization\$_{t-1}\$" D.wedis "$\Delta\$Party\$^\dag$ Polarization" L.wedis "Party\$^\dag$ Polarization\$_{t-1}\$" _cons Constant) substitute(\_ _) fragment
foreach i in 2 4{
esttab m`i'_wedis m`i'_polv m`i'_pi m`i'_estray m`i'_sartori using "Graphs/Appendix C`i'.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM se1_LRM b5_LRM se5_LRM b9_LRM se9_LRM, fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress ///
label equations(1,2) ///
rename(L.polv L.wedis L.pi L.wedis L.estray L.wedis L.sartori L.wedis D.polv D.wedis D.pi D.wedis D.estray D.wedis D.sartori D.wedis) ///
coeflabels(D.resleftvot "$\Delta\$Citizen Polarization (Voters)" L.resleftvot "Citizen Polarization\$_{t-1}\$ (Voters)" D.resleftabs "$\Delta\$Citizen Polarization (Non-voters)" L.resleftabs "Citizen Polarization\$_{t-1}\$ (Non-voters)" D.wedis "$\Delta\$Party\$^\dag$ Polarization" L.wedis "Party\$^\dag$ Polarization\$_{t-1}\$" _cons Constant) substitute(\_ _) fragment
}
foreach i in 3 5{
esttab m`i'_wedis m`i'_polv m`i'_pi m`i'_estray m`i'_sartori using "Graphs/Appendix C`i'.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM se1_LRM b5_LRM se5_LRM b9_LRM se9_LRM, fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress ///
label equations(1,2,3) ///
rename(L.polv L.wedis L.pi L.wedis L.estray L.wedis L.sartori L.wedis D.polv D.wedis D.pi D.wedis D.estray D.wedis D.sartori D.wedis) ///
coeflabels(D.resleftsop1 "$\Delta\$Citizen Polarization (High Sophistication)" L.resleftsop1 "Citizen Polarization\$_{t-1}\$ (High Sophistication)" D.resleftsop2 "$\Delta\$Citizen Polarization (Med. Sophistication)" L.resleftsop2 "Citizen Polarization\$_{t-1}\$ (Med. Sophistication)" ///
D.resleftsop3 "$\Delta\$Citizen Polarization (Low Sophistication)" L.resleftsop3 "Citizen Polarization\$_{t-1}\$ (Low Sophistication)" D.wedis "$\Delta\$Party\$^\dag$ Polarization" L.wedis "Party\$^\dag$ Polarization\$_{t-1}\$" _cons Constant) substitute(\_ _) fragment
}
*/

** Appendix D. SUR - Knowledgeable and Unknowledgeable Respondents
eststo d_m1: sureg (d.uwedis l.uwedis d.resleftpkn l.resleftpkn) (d.uwedis l.uwedis d.resleftpuk l.resleftpuk), dfk corr
eststo d_m2: sureg (d.resleftpkn l.resleftpkn d.uwedis l.uwedis) (d.resleftpuk l.resleftpuk d.uwedis l.uwedis), dfk corr

forval i=1/2{
est restore d_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 5{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
/*
esttab d_m1 d_m2 using "Graphs/Appendix D.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b5_LRM se1_LRM se5_LRM, layout(@ @ @ "@ @" "@ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.D1" "Model.D2") ///
order(L.uwedis D.resleftpkn D.resleftpuk L.resleftpkn L.resleftpuk D.uwedis) coeflabels(D.resleftpkn "$\Delta\$Citizen Polarization (Knowledgeable)" L.resleftpkn "Citizen Polarization\$_{t-1}\$ (Knowledgeable)" D.resleftpuk "$\Delta\$Citizen Polarization (Unknowledgeable)" L.resleftpuk "Citizen Polarization\$_{t-1}\$ (Unknowledgeable)" ///
D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" _cons Constant) substitute(\_ _) fragment
*/

** Appendix E. SUR - Highly and Lowly Educated Respondents
eststo e_m1: sureg (d.uwedis l.uwedis d.reslefteduch l.reslefteduch) (d.uwedis l.uwedis d.reslefteducl l.reslefteducl), dfk corr
eststo e_m2: sureg (d.reslefteduch l.reslefteduch d.uwedis l.uwedis) (d.reslefteducl l.reslefteducl d.uwedis l.uwedis), dfk corr

forval i=1/2{
est restore e_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 5{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
/*
esttab e_m1 e_m2 using "Graphs/Appendix E.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b5_LRM se1_LRM se5_LRM, layout(@ @ @ "@ @" "@ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.E1" "Model.E2") ///
order(L.uwedis D.reslefteduch D.reslefteducl L.reslefteduch L.reslefteducl D.uwedis) coeflabels(D.reslefteduch "$\Delta\$Citizen Polarization (High Education)" L.reslefteduch "Citizen Polarization\$_{t-1}\$ (High Education)" D.reslefteducl "$\Delta\$Citizen Polarization (Low Education)" L.reslefteducl "Citizen Polarization\$_{t-1}\$ (Low Education)" ///
D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" _cons Constant) substitute(\_ _) fragment
*/

** Appendix F1. SUR - wo/ gaps
preserve
keep if gap!=1
eststo f1_m1: sureg (d.uwedis l.uwedis d.resleft l.resleft) (d.resleft l.resleft d.uwedis l.uwedis), dfk corr
eststo f1_m2: sureg (d.uwedis l.uwedis d.resleftvot l.resleftvot) (d.uwedis l.uwedis d.resleftabs l.resleftabs), dfk corr
eststo f1_m3: sureg (d.uwedis l.uwedis d.resleftsop1 l.resleftsop1) (d.uwedis l.uwedis d.resleftsop2 l.resleftsop2) (d.uwedis l.uwedis d.resleftsop3 l.resleftsop3), dfk corr
eststo f1_m4: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis) (d.resleftabs l.resleftabs d.uwedis l.uwedis), dfk corr
eststo f1_m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis), dfk corr

forval i=1/5{
est restore f1_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 5 9{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
restore
/*
esttab f1_m1 f1_m2 f1_m3 f1_m4 f1_m5 using "Graphs/Appendix F1.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b5_LRM b9_LRM se1_LRM se5_LRM se9_LRM, layout(@ @ @ "@ @ @" "@ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) ///
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.F1.1" "Model.F1.2" "Model.F1.3" "Model.F1.4" "Model.F1.5") order(L.uwedis D.resleft L.resleft D.uwedis) /// 
rename(D.resleftvot D.resleft D.resleftabs D.resleft D.resleftsop1 D.resleft D.resleftsop2 D.resleft D.resleftsop3 D.resleft L.resleftvot L.resleft L.resleftabs L.resleft L.resleftsop1 L.resleft L.resleftsop2 L.resleft L.resleftsop3 L.resleft) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" _cons Constant) substitute(\_ _) fragment
*/

** Appendix F2. SUR - w/ time between two elections
eststo f2_m1: sureg (d.uwedis l.uwedis d.resleft l.resleft mdif) (d.resleft l.resleft d.uwedis l.uwedis mdif), dfk corr
eststo f2_m2: sureg (d.uwedis l.uwedis d.resleftvot l.resleftvot mdif) (d.uwedis l.uwedis d.resleftabs l.resleftabs mdif), dfk corr
eststo f2_m3: sureg (d.uwedis l.uwedis d.resleftsop1 l.resleftsop1 mdif) (d.uwedis l.uwedis d.resleftsop2 l.resleftsop2 mdif) (d.uwedis l.uwedis d.resleftsop3 l.resleftsop3 mdif), dfk corr
eststo f2_m4: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis mdif) (d.resleftabs l.resleftabs d.uwedis l.uwedis mdif), dfk corr
eststo f2_m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis mdif) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis mdif) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis mdif), dfk corr

forval i=1/5{
est restore f2_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 6 10{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
*disp "b`i'_LRM=" %9.3f b`i'_LRM
*disp "se`i'_LRM=" %9.3f se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
/*
esttab f2_m1 f2_m2 f2_m3 f2_m4 f2_m5 using "Graphs/Appendix F2.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b6_LRM b10_LRM se1_LRM se6_LRM se10_LRM, layout(@ @ @ "@ @ @" "@ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) ///
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.F2.1" "Model.F2.2" "Model.F2.3" "Model.F2.4" "Model.F2.5") order(L.uwedis D.resleft L.resleft D.uwedis) /// 
rename(D.resleftvot D.resleft D.resleftabs D.resleft D.resleftsop1 D.resleft D.resleftsop2 D.resleft D.resleftsop3 D.resleft L.resleftvot L.resleft L.resleftabs L.resleft L.resleftsop1 L.resleft L.resleftsop2 L.resleft L.resleftsop3 L.resleft) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" mdif "Time bw. Elections (in Months)" _cons Constant) substitute(\_ _) fragment
*/

** Appendix G. SUR - w/ controls for election decade
eststo g_m1: sureg (d.uwedis l.uwedis d.resleft l.resleft i.decade) (d.resleft l.resleft d.uwedis l.uwedis i.decade), dfk corr
eststo g_m2: sureg (d.uwedis l.uwedis d.resleftvot l.resleftvot i.decade) (d.uwedis l.uwedis d.resleftabs l.resleftabs i.decade), dfk corr
eststo g_m3: sureg (d.uwedis l.uwedis d.resleftsop1 l.resleftsop1 i.decade) (d.uwedis l.uwedis d.resleftsop2 l.resleftsop2 i.decade) (d.uwedis l.uwedis d.resleftsop3 l.resleftsop3 i.decade), dfk corr
eststo g_m4: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis i.decade) (d.resleftabs l.resleftabs d.uwedis l.uwedis i.decade), dfk corr
eststo g_m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis i.decade) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis i.decade) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis i.decade), dfk corr

forval i=1/5{
est restore g_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 10 19{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
/*
esttab g_m1 g_m2 g_m3 g_m4 g_m5 using "Graphs/Appendix G.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b10_LRM b19_LRM se1_LRM se10_LRM se19_LRM, layout(@ @ @ "@ @ @" "@ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) /// 
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.G1" "Model.G2" "Model.G3" "Model.G4" "Model.G5") order(L.uwedis D.resleft L.resleft D.uwedis) /// 
rename(D.resleftvot D.resleft D.resleftabs D.resleft D.resleftsop1 D.resleft D.resleftsop2 D.resleft D.resleftsop3 D.resleft L.resleftvot L.resleft L.resleftabs L.resleft L.resleftsop1 L.resleft L.resleftsop2 L.resleft L.resleftsop3 L.resleft) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" 1.decade "Election Decade=1980s" 2.decade "Election Decade=1990s" ///
3.decade "Election Decade=2000s" 4.decade "Election Decade=2010s" _cons Constant) substitute(\_ _) fragment
*/

** Appendix H. SUR - w/ ENEP
eststo h_m1: sureg (d.uwedis l.uwedis d.resleft l.resleft d.enep l.enep) (d.resleft l.resleft d.uwedis l.uwedis d.enep l.enep), dfk corr
eststo h_m2: sureg (d.uwedis l.uwedis d.resleftvot l.resleftvot d.enep l.enep) (d.uwedis l.uwedis d.resleftabs l.resleftabs d.enep l.enep), dfk corr
eststo h_m3: sureg (d.uwedis l.uwedis d.resleftsop1 l.resleftsop1 d.enep l.enep) (d.uwedis l.uwedis d.resleftsop2 l.resleftsop2 d.enep l.enep) (d.uwedis l.uwedis d.resleftsop3 l.resleftsop3 d.enep l.enep), dfk corr
eststo h_m4: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis d.enep l.enep) (d.resleftabs l.resleftabs d.uwedis l.uwedis d.enep l.enep), dfk corr
eststo h_m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis d.enep l.enep) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis d.enep l.enep) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis d.enep l.enep), dfk corr

forval i=1/5{
est restore h_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 7 13{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
/*
esttab h_m1 h_m2 h_m3 h_m4 h_m5 using "Graphs/Appendix H.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b7_LRM b13_LRM se1_LRM se7_LRM se13_LRM, layout(@ @ @ "@ @ @" "@ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) /// 
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.H1" "Model.H2" "Model.H3" "Model.H4" "Model.H5") order(L.uwedis D.resleft L.resleft D.uwedis D.enep L.enep) /// 
rename(D.resleftvot D.resleft D.resleftabs D.resleft D.resleftsop1 D.resleft D.resleftsop2 D.resleft D.resleftsop3 D.resleft L.resleftvot L.resleft L.resleftabs L.resleft L.resleftsop1 L.resleft L.resleftsop2 L.resleft L.resleftsop3 L.resleft) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" L.enep "Effective Number of Parties\$_{t-1}\$" ///
D.enep "$\Delta\$Effective Number of Parties" _cons Constant) substitute(\_ _) fragment
*/

** Appendix I. SUR - w/ average district magnitude
eststo i_m1: sureg (d.uwedis l.uwedis d.resleft l.resleft d.av_dismag l.av_dismag) (d.resleft l.resleft d.uwedis l.uwedis d.av_dismag l.av_dismag), dfk corr
eststo i_m2: sureg (d.uwedis l.uwedis d.resleftvot l.resleftvot d.av_dismag l.av_dismag) (d.uwedis l.uwedis d.resleftabs l.resleftabs d.av_dismag l.av_dismag), dfk corr
eststo i_m3: sureg (d.uwedis l.uwedis d.resleftsop1 l.resleftsop1 d.av_dismag l.av_dismag) (d.uwedis l.uwedis d.resleftsop2 l.resleftsop2 d.av_dismag l.av_dismag) (d.uwedis l.uwedis d.resleftsop3 l.resleftsop3 d.av_dismag l.av_dismag), dfk corr
eststo i_m4: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis d.av_dismag l.av_dismag) (d.resleftabs l.resleftabs d.uwedis l.uwedis d.av_dismag l.av_dismag), dfk corr
eststo i_m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis d.av_dismag l.av_dismag) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis d.av_dismag l.av_dismag) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis d.av_dismag l.av_dismag), dfk corr

forval i=1/5{
est restore i_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 7 13{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
/*
esttab i_m1 i_m2 i_m3 i_m4 i_m5 using "Graphs/Appendix I.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b7_LRM b13_LRM se1_LRM se7_LRM se13_LRM, layout(@ @ @ "@ @ @" "@ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) ///
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.I1" "Model.I2" "Model.I3" "Model.I4" "Model.I5") order(L.uwedis D.resleft L.resleft D.uwedis D.av_dismag L.av_dismag) ///  
rename(D.resleftvot D.resleft D.resleftabs D.resleft D.resleftsop1 D.resleft D.resleftsop2 D.resleft D.resleftsop3 D.resleft L.resleftvot L.resleft L.resleftabs L.resleft L.resleftsop1 L.resleft L.resleftsop2 L.resleft L.resleftsop3 L.resleft) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" L.av_dismag "Average District Magnitude\$_{t-1}\$" ///
D.av_dismag "$\Delta\$Average District Magnitude" _cons Constant) substitute(\_ _) fragment
*/

** Appendix J1. SUR - w/ controls for electoral rule
eststo j1_m1: sureg (d.uwedis l.uwedis d.resleft l.resleft i.pr_rule) (d.resleft l.resleft d.uwedis l.uwedis i.pr_rule), dfk corr
eststo j1_m2: sureg (d.uwedis l.uwedis d.resleftvot l.resleftvot i.pr_rule) (d.uwedis l.uwedis d.resleftabs l.resleftabs i.pr_rule), dfk corr
eststo j1_m3: sureg (d.uwedis l.uwedis d.resleftsop1 l.resleftsop1 i.pr_rule) (d.uwedis l.uwedis d.resleftsop2 l.resleftsop2 i.pr_rule) (d.uwedis l.uwedis d.resleftsop3 l.resleftsop3 i.pr_rule), dfk corr
eststo j1_m4: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis i.pr_rule) (d.resleftabs l.resleftabs d.uwedis l.uwedis i.pr_rule), dfk corr
eststo j1_m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis i.pr_rule) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis i.pr_rule) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis i.pr_rule), dfk corr

forval i=1/5{
est restore j1_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 8 15{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
/*
esttab j1_m1 j1_m2 j1_m3 j1_m4 j1_m5 using "Graphs/Appendix J1.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b8_LRM b15_LRM se1_LRM se8_LRM se15_LRM, layout(@ @ @ "@ @ @" "@ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) ///
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.J1.1" "Model.J1.2" "Model.J1.3" "Model.J1.4" "Model.J1.5") order(L.uwedis D.resleft L.resleft D.uwedis) ///  
rename(D.resleftvot D.resleft D.resleftabs D.resleft D.resleftsop1 D.resleft D.resleftsop2 D.resleft D.resleftsop3 D.resleft L.resleftvot L.resleft L.resleftabs L.resleft L.resleftsop1 L.resleft L.resleftsop2 L.resleft L.resleftsop3 L.resleft) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" 1.pr_rule "Mixed" 2.pr_rule "Majoritarian" _cons Constant) substitute(\_ _) fragment
*/

** Appendix J2. SUR - SMD interaction
gen smd=pr_rule==2
eststo j2_m1: sureg (d.uwedis l.uwedis d.resleft l.resleft smd c.l.uwedis#smd c.d.resleft#smd c.l.resleft#smd) (d.resleft l.resleft d.uwedis l.uwedis smd c.l.resleft#smd c.d.uwedis#smd c.l.uwedis#smd), dfk corr
eststo j2_m2: sureg (d.uwedis l.uwedis d.resleftvot l.resleftvot smd c.l.uwedis#smd c.d.resleftvot#smd c.l.resleftvot#smd) (d.uwedis l.uwedis d.resleftabs l.resleftabs smd c.l.uwedis#smd c.d.resleftabs#smd c.l.resleftabs#smd), dfk corr
eststo j2_m3: sureg (d.uwedis l.uwedis d.resleftsop1 l.resleftsop1 smd c.l.uwedis#smd c.d.resleftsop1#smd c.l.resleftsop1#smd) (d.uwedis l.uwedis d.resleftsop2 l.resleftsop2 smd c.l.uwedis#smd c.d.resleftsop2#smd c.l.resleftsop2#smd) (d.uwedis l.uwedis d.resleftsop3 l.resleftsop3 smd c.l.uwedis#smd c.d.resleftsop3#smd c.l.resleftsop3#smd), dfk corr
eststo j2_m4: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis smd c.l.resleftvot#smd c.d.uwedis#smd c.l.uwedis#smd) (d.resleftabs l.resleftabs d.uwedis l.uwedis smd c.l.resleftabs#smd c.d.uwedis#smd c.l.uwedis#smd), dfk corr
eststo j2_m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis smd c.l.resleftsop1#smd c.d.uwedis#smd c.l.uwedis#smd) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis smd c.l.resleftsop2#smd c.d.uwedis#smd c.l.uwedis#smd) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis smd c.l.resleftsop3#smd c.d.uwedis#smd c.l.uwedis#smd), dfk corr

/*
esttab j2_m1 j2_m2 j2_m3 j2_m4 j2_m5 using "Graphs/Appendix J2.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b5_LRM b9_LRM se1_LRM se5_LRM se9_LRM, layout(@ @ @ "@ @ @" "@ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) ///
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.J2.1" "Model.J2.2" "Model.J2.3" "Model.J2.4" "Model.J2.5") order(L.uwedis D.resleft L.resleft D.uwedis) ///
rename(D.resleftvot D.resleft D.resleftabs D.resleft D.resleftsop1 D.resleft D.resleftsop2 D.resleft D.resleftsop3 D.resleft L.resleftvot L.resleft L.resleftabs L.resleft L.resleftsop1 L.resleft L.resleftsop2 L.resleft L.resleftsop3 L.resleft  ///
1.smd#cD.resleftvot 1.smd#D.resleft 1.smd#cD.resleftabs 1.smd#D.resleft 1.smd#cD.resleftsop1 1.smd#D.resleft 1.smd#cD.resleftsop2 1.smd#D.resleft 1.smd#cD.resleftsop3 1.smd#D.resleft 1.smd#cL.resleftvot 1.smd#L.resleft 1.smd#cL.resleftabs ///
1.smd#L.resleft 1.smd#cL.resleftsop1 1.smd#L.resleft 1.smd#cL.resleftsop2 1.smd#L.resleft 1.smd#cL.resleftsop3 1.smd#L.resleft) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" ///
1.smd#cL.resleft "Majoritarian $\times$ Citizen\$^\dag$ Polarization\$_{t-1}\$" 1.smd#cD.resleft "Majoritarian $\times$ $\Delta$ Citizen\$^\dag$ Polarization" 1.smd#cL.uwedis "Majoritarian $\times$ Party Polarization\$_{t-1}\$" ///
1.smd#cD.uwedis "Majoritarian $\times$ $\Delta$ Party Polarization" smd "Majoritarian" _cons Constant) substitute(\_ _) fragment
*/

** Appendix K - OLS w/ clustered standard errors by party systems, Prais-Winsten, and Random-effects GLM regressions
eststo ols: reg d.uwedis l.uwedis d.resleft l.resleft, robust cluster(iso)
eststo olsb: reg d.resleft l.resleft d.uwedis l.uwedis, robust cluster(iso)
eststo xtpcse: xtpcse d.uwedis l.uwedis d.resleft l.resleft, correlation(ar1) np1 rhotype(tscorr) pairwise
eststo xtpcseb: xtpcse d.resleft l.resleft d.uwedis l.uwedis, correlation(ar1) np1 rhotype(tscorr) pairwise
eststo xtregre: xtreg d.uwedis l.uwedis d.resleft l.resleft, re
eststo xtregreb: xtreg d.resleft l.resleft d.uwedis l.uwedis, re

foreach m in ols olsb xtpcse xtpcseb xtregre xtregreb {
est restore `m'
mat b=e(b)
mat V=e(V)
sca df=e(N)-e(df_m)-1
loc i=1
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(df, t`i'_LRM)
}
/*
esttab ols olsb xtpcse xtpcseb xtregre xtregreb using "Graphs/Appendix K.tex", tex replace b(%10.3f) se stats(N r2 r2_o b1_LRM se1_LRM, fmt(0 3) labels("N" "\$R^{2}\$" "\$R^2_{overall}\$")) ///
label starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.K1" "Model.K2" "Model.K3" "Model.K4" "Model.K5" "Model.K6" "Model.K7" "Model.K8") order(L.uwedis D.resleft L.resleft D.uwedis) ///
coeflabels(D.resleft "$\Delta\$Citizen Polarization (All Respondents)" L.resleft "Citizen Polarization\$_{t-1}\$ (All Respondents)" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" _cons Constant) substitute(\_ _) fragment
*/

** Appendix L - SUR w/ 2-election lagged DVs
eststo l_m1: sureg (d.uwedis l.uwedis d.resleft l.resleft l2.uwedis) (d.resleft l.resleft d.uwedis l.uwedis l2.resleft), dfk corr
eststo l_m2: sureg (d.uwedis l.uwedis d.resleftvot l.resleftvot l2.uwedis) (d.uwedis l.uwedis d.resleftabs l.resleftabs l2.uwedis), dfk corr
eststo l_m3: sureg (d.uwedis l.uwedis d.resleftsop1 l.resleftsop1 l2.uwedis) (d.uwedis l.uwedis d.resleftsop2 l.resleftsop2 l2.uwedis) (d.uwedis l.uwedis d.resleftsop3 l.resleftsop3 l2.uwedis), dfk corr
eststo l_m4: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis l2.resleftvot) (d.resleftabs l.resleftabs d.uwedis l.uwedis l2.resleftabs), dfk corr
eststo l_m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis l2.resleftsop1) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis l2.resleftsop2) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis l2.resleftsop3), dfk corr

forval i=1/5{
est restore l_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 6 10{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
/*
esttab l_m1 l_m2 l_m3 l_m4 l_m5 using "Graphs/Appendix L.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b6_LRM b10_LRM se1_LRM se6_LRM se10_LRM, layout(@ @ @ "@ @ @" "@ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) ///
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.L1" "Model.L2" "Model.L3" "Model.L4" "Model.L5") order(L.uwedis D.resleft L.resleft D.uwedis) ///
rename(D.resleftvot D.resleft D.resleftabs D.resleft D.resleftsop1 D.resleft D.resleftsop2 D.resleft D.resleftsop3 D.resleft L.resleftvot L.resleft L.resleftabs L.resleft L.resleftsop1 L.resleft L.resleftsop2 L.resleft L.resleftsop3 L.resleft ///
L2.resleft L2.uwedis L2.resleftvot L2.uwedis L2.resleftabs L2.uwedis L2.resleftsop1 L2.uwedis L2.resleftsop2 L2.uwedis L2.resleftsop3 L2.uwedis) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" L2.uwedis "Dependent Variable\$_{t-2}\$" _cons Constant) substitute(\_ _) fragment
*/

** Appendix M - Other Scales
eststo m_m1: sureg (d.uwedis l.uwedis d.resleft l.resleft i.other_scale) (d.resleft l.resleft d.uwedis l.uwedis i.other_scale), dfk corr
eststo m_m2: sureg (d.uwedis l.uwedis d.resleftvot l.resleftvot i.other_scale) (d.uwedis l.uwedis d.resleftabs l.resleftabs i.other_scale), dfk corr
eststo m_m3: sureg (d.uwedis l.uwedis d.resleftsop1 l.resleftsop1 i.other_scale) (d.uwedis l.uwedis d.resleftsop2 l.resleftsop2 i.other_scale) (d.uwedis l.uwedis d.resleftsop3 l.resleftsop3 i.other_scale), dfk corr
eststo m_m4: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis i.other_scale) (d.resleftabs l.resleftabs d.uwedis l.uwedis i.other_scale), dfk corr
eststo m_m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis i.other_scale) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis i.other_scale) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis i.other_scale), dfk corr

forval i=1/5{
est restore m_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 11 21{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}

/*
esttab m_m1 m_m2 m_m3 m_m4 m_m5 using "Graphs/Appendix M.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b11_LRM b21_LRM se1_LRM se11_LRM se21_LRM, layout(@ @ @ "@ @ @" "@ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) ///
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.C1" "Model.C2" "Model.C3" "Model.C4" "Model.C5") order(L.uwedis D.resleft L.resleft D.uwedis) ///
rename(D.resleftvot D.resleft D.resleftabs D.resleft D.resleftsop1 D.resleft D.resleftsop2 D.resleft D.resleftsop3 D.resleft L.resleftvot L.resleft L.resleftabs L.resleft L.resleftsop1 L.resleft L.resleftsop2 L.resleft L.resleftsop3 L.resleft) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" 2.other_scale "0/6 Scale" 3.other_scale "1/10 Scale" ///
4.other_scale "1/11 Scale" 5.other_scale "1/7 Scale" 6.other_scale "1/9 Scale" _cons Constant) substitute(\_ _) fragment
*/
** Appendix N - SUR w/ compulsory voting
eststo n_m1: sureg (d.uwedis l.uwedis d.resleft l.resleft compulsory) (d.resleft l.resleft d.uwedis l.uwedis compulsory), dfk corr
eststo n_m2: sureg (d.uwedis l.uwedis d.resleftvot l.resleftvot compulsory) (d.uwedis l.uwedis d.resleftabs l.resleftabs compulsory), dfk corr
eststo n_m3: sureg (d.uwedis l.uwedis d.resleftsop1 l.resleftsop1 compulsory) (d.uwedis l.uwedis d.resleftsop2 l.resleftsop2 compulsory) (d.uwedis l.uwedis d.resleftsop3 l.resleftsop3 compulsory), dfk corr
eststo n_m4: sureg (d.resleftvot l.resleftvot d.uwedis l.uwedis compulsory) (d.resleftabs l.resleftabs d.uwedis l.uwedis compulsory), dfk corr
eststo n_m5: sureg (d.resleftsop1 l.resleftsop1 d.uwedis l.uwedis compulsory) (d.resleftsop2 l.resleftsop2 d.uwedis l.uwedis compulsory) (d.resleftsop3 l.resleftsop3 d.uwedis l.uwedis compulsory), dfk corr

forval i=1/5{
est restore n_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 6 10{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
/*
esttab n_m1 n_m2 n_m3 n_m4 n_m5 using "Graphs/Appendix N.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b6_LRM b10_LRM se1_LRM se6_LRM se10_LRM, layout(@ @ @ "@ @ @" "@ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) ///
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.N1" "Model.N2" "Model.N3" "Model.N4" "Model.N5") order(L.uwedis D.resleft L.resleft D.uwedis) ///
rename(D.resleftvot D.resleft D.resleftabs D.resleft D.resleftsop1 D.resleft D.resleftsop2 D.resleft D.resleftsop3 D.resleft L.resleftvot L.resleft L.resleftabs L.resleft L.resleftsop1 L.resleft L.resleftsop2 L.resleft L.resleftsop3 L.resleft) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" compulsory "Compulsory Voting" _cons Constant) substitute(\_ _) fragment
*/

** Appendix O. SUR - Highest, Med. High, Med.Low, and Lowest Sophistication
eststo o_m1: sureg (d.uwedis l.uwedis d.resleftsop_alt1 l.resleftsop_alt1) (d.uwedis l.uwedis d.resleftsop_alt2 l.resleftsop_alt2) (d.uwedis l.uwedis d.resleftsop_alt3 l.resleftsop_alt3) (d.uwedis l.uwedis d.resleftsop_alt4 l.resleftsop_alt4), dfk corr
eststo o_m2: sureg (d.resleftsop_alt1 l.resleftsop_alt1 d.uwedis l.uwedis) (d.resleftsop_alt2 l.resleftsop_alt2 d.uwedis l.uwedis) (d.resleftsop_alt3 l.resleftsop_alt3 d.uwedis l.uwedis) (d.resleftsop_alt4 l.resleftsop_alt4 d.uwedis l.uwedis), dfk corr

forval i=1/2{
est restore o_m`i'
mat b=e(b)
mat V=e(V)
foreach i of numlist 1 5 9 13{
sca x`i'=b[1,`=`i'+2']
sca y`i'=b[1,`i']
sca varx`i'=V[`=`i'+2',`=`i'+2']
sca vary`i'=V[`i',`i']
sca covx`i'y`i'=V[`=`i'+2',`i']
sca b`i'_LRM=x`i'/-y`i'
sca se`i'_LRM=sqrt((1/y`i'^2)*varx`i'+(x`i'^2/y`i'^4)*vary`i'-2*(x`i'/y`i'^3)*covx`i'y`i')
sca t`i'_LRM=b`i'_LRM/se`i'_LRM
estadd sca b`i'_LRM
estadd sca se`i'_LRM
disp "p`i'_LRM=" %9.3f ttail(e(dfk2_adj), t`i'_LRM)
}
}
/*
esttab o_m1 o_m2 using "Graphs/Appendix O.tex", tex replace b(%10.3f) se stats(N r2 chi2_bp b1_LRM b5_LRM b9_LRM b13_LRM se1_LRM se5_LRM se9_LRM se13_LRM, layout(@ @ @ "@ @ @ @" "@ @ @ @") fmt(0 3 3) labels("N" "\$R^{2}\$" "Breusch-Pagan Test\$^\ddag\$ (\$\chi^{2}\$)")) ///
unstack starlevels(* 0.05 ** 0.01) alignment(l) nonote nogaps nonumbers noomit nobase compress mtitle("Model.O1" "Model.O2") order(L.uwedis D.resleft L.resleft D.uwedis) /// 
rename(D.resleftsop_alt1 D.resleft D.resleftsop_alt2 D.resleft D.resleftsop_alt3 D.resleft D.resleftsop_alt4 D.resleft L.resleftsop_alt1 L.resleft L.resleftsop_alt2 L.resleft L.resleftsop_alt3 L.resleft L.resleftsop_alt4 L.resleft) ///
coeflabels(D.resleft "$\Delta\$Citizen\$^\dag$ Polarization" L.resleft "Citizen\$^\dag$ Polarization\$_{t-1}\$" D.uwedis "$\Delta\$Party Polarization" L.uwedis "Party Polarization\$_{t-1}\$" _cons Constant) substitute(\_ _) fragment
*/
